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We study the topological susceptibility, \, in two flavour lattice QCD.EJ We find 
OO ■ clear evidence for the expected suppression of \ at small quark mass. The estimate 

of the pion decay constant, /„- = 105 ± 5 _,„ MeV, is consistent with the experi- 
mental value of approximately 93 MeV. We compare \ to the large- N c prediction 
^- ' and find consistency over a large range of quark masses. 

OO ' 

o 

^^ ' The ability to access the non-perturbative sectors, and to vary parameters 

0^ , fixed in Nature has made lattice Monte Carlo simulation a valuable tool fi 

o 

f— ») , The topological susceptibility is the squared expectation value of the topo 

logical charge, normalised by the volume 



- * — i 

X 



(Q 2 ) „_ i f*J. 
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X = ^-, <3= 32^7 ^x-e^ T F^(x)F: T (x)- (1) 

(-h _ Sea quarks induce an instanton-anti-instanton attraction which in the chiral 

limit becomes stronger, suppressing Q and \" 



X — S (m u + rrid ) , where £ = — lim lim (Q\yjij}\Q) (2) 

is the chiral condensate, a We assume (0|V>?/>|0) = (0|uw|0) = (0\dd\0) and ne- 
glect contributions of heavier quarks. The Gell-Mann-Oakes-Renner relation, 



tfml = (m u +m d )X+0(m 2 q ) =► x = J -^+0{m\) (3) 






for Nf light flavours, in a convention where the experimental value of the pion 
decay constant / ff ~ 93 MeV. Eq. || holds in the limit f^m^V » 1, which is 
satisfied by all our lattices. The higher order terms ensure that \ ~~> X q "\ the 
quenched value, as m q , m, — ► oo. We find, however, that our measured values 
are not very much smaller than x qu , so we must consider two possibilities. 
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Figure 1: The measured topological susceptibility, with interpolated quenched points at 
the same rrj and fits independent of the quenched points. 



-- fit to Eqn. 6 

-- fit to Eqn. 7 
— fit to Eqn. 
-- fit to Eqn. 9 
* physical value 
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Firstly, there are phenomenological reasons Lm for believing that QCD is 
'close' to N c = co_and in the case of gluodynamics even SU(2) is demonstrably 



q u 



close to SU(oo). Qfl Fermion effects are non-leading in N c , so we expect \ — ► x 
for any fixed value of m q as the number of colours iV c — > oo. For small 



/oo the quantities at leading order 
and perhaps higher order terms are 



m q we-jexpect Eq. 4 to hold, with x°° 

in N c . u Alternatively, our m q ~ ™ s trang 

important. In the absence of a QCD prediction, Eq. 5 interpolates between 

Eq. and the quenched limit. 
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— m^arctanf -J~x q " 
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(5). 



Measurements of x were made on a number of ensembles of Nf — 2 lattice 
field configurations produced by the UKQCD Collaboration n An SU(3) Wilson 
gauge action is coupled to clover improved Wilson fermions. u The UKQCD en- 
sembles have two notable features. The improvement is fully non-perturbative, 
with discretisation errors being quadratic rather than linear in the lattice spac- 
ing. Second, the couplings are chosen to maintain an approximately constant 
lattice spacing (as defined by the Sommer scale, tq — 0.49 fmEl) as the quark 
mass is varied. This is important, as the susceptibility in gluodynamics varies 
with the lattice spacing asH f^x — 0.072 — 0.208/rQ in competition with the 
variation with m q . u The topological susceptibility is measured from the gauge 
fields after cooling to remove the UV noise. We plot these data in Fig. [j] along 
with the interpolated x qu at an equivalent lattice spacing from the above for- 
mula for comparison, which vary little owing to the the UKQCD matching. The 
behaviour with M = (rom^) 2 is qualitatively as expected and, more quantita- 



Table 1: Fits to the JVfj t most chiral points of \. 



Fit 


N Rt 


CO 


c\ 


X 2 /d.o.f. 


fofn 


Eq. 6 
Eq. 6 
Eq. 6 


2 
3 

4 


0.0140 (16) 
0.0112 (6) 
0.0091 (4) 


— 


0.805 
2.202 
9.008 


0.237 (14) 
0.212 (6) 


Eq. 7 
Eq. 7 
Eq. 7 


3 

4 
5 


0.0176 (35) (4) 
0.0170 (16) (1) 
0.0147 (14) (1) 


-0.0018 (10) (1) 
-0.0016 (4) (0) 
-0.0011 (3) (0) 


0.964 
0.502 
2.965 


0.265 (27) 
0.261 (13) 

0.242 (12) 


Fit 


N Rt 


CO 


C3 


X 2 /d.o.f. 


fofn 


Eq. 8 
Eq. 8 
Eq. 8 


3 

4 
5 


0.0208 (87) (12) 
0.0272 (85) (18) 
0.0233 (66) (10) 


0.0844 (427) (35) 
0.0632 (114) (6) 
0.0717 (147) (3) 


1.013 
0.895 

1.847 


0.288 (61) 
0.329 (53) 
0.305 (44) 


Eq. 9 
Eq. 9 
Eq. 9 


3 

4 
5 


0.0186 (53) (7) 
0.0209 (42) (7) 
0.0189 (36) (5) 


0.0576 (175) (6) 
0.0506 (55) (5) 
0.0550 (69) (6) 


0.990 
0.682 
1.929 


0.273 (40) 
0.289 (30) 
0.275 (27) 



tively, we attempt fits motivated by Eqs. [|, 4, 5 

C0C3 



rp 2 X 
M 



co (6), co + ci(rom,r) (7), 



C3+C M W ' 7T \2C MJ 



We include progressively less chiral points until the fit becomes unacceptably 
bad in Table yj. We note the wide range fitted simply by including an m\ term, 
and the consistency of our data with large- N c predictions. The stability and 
similarity of the fits motivates us to use Co from Eq. 7 to estimate f T = 105 ± 
5 _iq MeV, with variation between other fits providing the second, systematic 
error, and in good agreement with the experimental value ~ 93 MeV. 
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